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ABSTRACT 

To calibrate stellar velocity dispersion measurements from optical and near-IR stellar lines, and to 
improve the black hole mass (Mbh)- stellar velocity dispersion (<j*) relation, we measure c* based on 
high quality iJ-band spectra for a sample of 31 nearby galaxies, for which dynamical A^bh is available 
in the literature. By comparing velocity dispersions measured from stellar lines in the -H-band with 
those measured from optical stellar lines, we find no significant difference, suggesting that optical and 
near-IR stellar lines represent the same kinematics and that dust effect is negligible for early-type 
galaxies. Based on the spatially-resolved rotation and velocity dispersion measurements along the 
major axis of each galaxy, we find that a rotating stellar disk is present for 80% of galaxies in the 
sample. For galaxies with a rotation component, <r* measured from a single aperture spectrum can 
vary by up to ~20%, depending on the size of the adopted extraction aperture. To correct for the 
rotational broadening, we derive luminosity-weighted er* within the effective radius of each galaxy, 
providing uniformly measured velocity dispersions to improve the Mbh-c* relation. 
Subject headings: galaxies: kinematics and dynamics-galaxies: bulges-infrared: galaxies-techniques: 
spectroscopic 
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1. INTRODUCTION 

The black hole mass (Mbh) correlation with host 
galaxy properties has been one of the main issues in un- 
derstanding galaxy evolution and black hole growth. In 
particular, the relatively tight correlation between A/bh 
and stellar velocity dispersion (Mbh-c*) has been re- 
porte d for nearby galaxies with dynamically measure d 
Mbh (jFerrarese fc Merritdl200"cl iGebhardt et al.ll2000ah . 
as well as present-day active galaxies with Mbh deter- 
mine d from reverberati on-mapping results (|Onken et al.l 
12004 IWoo et al.l I2O1O0 . While early studies claimed 
a remarkably tight Mbh-c* relation with its int rinsic 
scatter below ^0.3 dex (e.g.. iTremaine et al.l 120021 ) . re- 
cent studies showed a larger intrinsic scatter and a 
steeper slope due to the increased sample size, inclu- 
sion of more diverse galaxies, i.e., late- type and pseudo 
bulge galaxies, and the improvements of Mbh measurc- 
ments based on better dynamical modeling and data 
dFerrarese fc Ford! [2005 iGrahaml [200l iGultekin et all 
2009at IMcConnell et al ] I2011LI2012D ~ 

In understanding BH-galaxy coevolution, the present- 
day Mbh-c* relation sets a local calibration point as most 
observational studies investigated cosmic evolution of the 
Mbh-c* relation by measuring an offset from the local 
relationship (e.g., IWoo et al.l [20061 120081 : iBennert et all 
I2010L 1201 lbf ). At the same time, the present-day M B h- 
cr* relation has been used for calibrating the Mbh of 
active galactic nuclei (AGN), which is determined from 
the kinematics of sub-pc scale broad-emission line region. 

7 Author to whom any correspondence should be addressed 



The unknown viral factor for converting the line-of-sight 
velocity of broad-line region gas to the intrinsic veloc- 
ity, has been empirically determined by matching the 
Mbh-c* relati on of quiescent an d acti ve galaxies a t z ~0 
(IQnken et al.l l200l IWolTetlll f2?Tint IPark et al.ll20Tl . 
Thus, defining the Mbh-c* relation in the local universe 
is of importance to unveil the nature of BH-galaxy co- 
evolution. 

Stellar kinematics studies based on the near-IR stel- 
lar lines became powerful as near-IR spectrographs com- 
bined with laser-guide star adaptive optics provides the 
best spatial resolution f or the ground-based facilities 
(e.g.. iWatson et all 120081 ). Moreover, measuring c* in 
the near-IR is more promising for AGN host galaxies 
since AGN-to-s tar flux ratios are much more favorable 
in the near-IR (iDasvra et al.l 120071 IWatson et all 120081 : 
IWoo et al.l I2010T ) while it is almost impossible to mea- 
sure a* in the optical for host galaxies of high luminosity 
QSOs. 

Despite the increasing usage of near-IR spectra for 
probing stellar kinematics, a proper comparison between 
optical and near-IR measurements is still lacking. By 
measuring velocity dispersion of 25 early-type galaxies 
based on the CO absorption band head at 2.29/xm in the 
iGband, iSilge fc Gebhardtl (|2003ft claimed that velocity 
dispersion measured from near-IR stellar lines was sys- 
tematically smaller by 10-30% than that measured from 
optica l stellar lines. In contrast, iRothberg fc Fischer! 
(2010) reported that optical and near-IR velocity dis- 
persions were consistent for a sample of 23 early-type 
galaxies, by comparing a* measured from the CO band 
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heads i n the if -band, with cr* m easured from the Call 
triplet. Vandcrbek e et all (|2011l ) also measured velocity 
dispersion based on the CO band heads for a sample of 
22 early-type galaxies, and presented consistent results 
with respect to optical cr*. The discrepancy among var- 
ious studies may have resulted from the systematic un- 
certainties of the velocity dispersion measurements since 
the line dispersion was measured from intrinsically broad 
CO band heads in the if- band and template mism atch 
could be very strong (see iSilge fc GebhardtJ 120031 ) . In 
contrast, the ii-band spectral range (~1.6-~1.7jUm) con- 
tains many more stellar lines, e.g., Si I, CO, and Mg I 
than the if-band, and is possibly less susceptible to tem- 
plate mismatch although the presence of strong sky OH 
lines is a downside. To utilize the ii-band stellar lines 
for studying stellar kinematics, a proper comparison is re- 
quired between cr* measured from ii-band spectra with 
that measured from optical spectra stellar lines. 

To derive reliable cr* to represent the kinematics of 
the pressure-supported bulge or spheroidal component, 
the effect of the rotation component should be corrected 
for. In the case of galaxies with a rotating stellar disk, 
the line-of-sight velocity dispersion can be easily overes- 
timated due to rotational broadeni ng if a large aperture 
is used to extract a spectrum (e.g.. iBennert et a"Lll2011at 
lHarris et al.ll2012l) . The effect of rotational broadening is 
stronger for more edge-on stellar disks, potentially pro- 
ducing systematic bias. Thus, it is important to correct 
for rotational broadening. 

In this paper, we measure the stellar velocity disper- 
sion of 31 nearby galaxies using high quality ii-band 
spectra. We compare cr* measurements based on stel- 
lar lines in the ii-band with optical cr* measurements 
from the literature. We also correct for the rotation and 
aperture effect based on the spatially resolved kinemat- 
ics measurements to improve the Mbh-c* relation. The 
paper is organized as follows. We describe sample selec- 
tion, observations and data reduction in § [21 In § 3, we 
present cr* measurements, and the effects of rotation and 
aperture size. In § 4, we compare our ii-band cr* mea- 
surements with optical cr* from the literature, and derive 
the Mbh-c* relation for early-type galaxies based on the 
rotation-corrected cr* . The main results are summarized 
in §[5] 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection and Observations 

To directly compare optical and near-IR stellar veloc- 
ity dispersions and to calibrate the Mbh-c* relation, we 
select 31 nearby galax ies from the Mbh-c* sample (e.g., 
iGultekin et al.l l2009af l . for which dynamical AiBH mea- 
surements and optical cr* measurements are available. 
The sample is mainly composed of early-type galaxies (20 
ellipticals, 8 lenticulars and 3 spirals) and spans a wide 
range in cr* from 67 km s^ 1 to 385 km s' 1 as listed in Ta- 
ble [T] Also, it covers three orders of magnitude in Mbh 
and constitutes about half the sam ple size of galaxies 
with dynamical Mbh measurements (jMcConnell fc Ma"l 
l20ll . 

Observations were performed at the Palomar Hale 5 m 
telescope using the near-IR spectrograph TripleSpec, si- 
multaneously covering the wavelength range from 1.0 /jm 
to 2.4 /im. In this work we only employ the ii-band spec- 



tra centered at ~1.7 /^m, as it covers many stellar absorp- 
tion lines suitable for the cr* measurement. We place an 
1" x 30" long-slit along the major axis of each galaxy. 
The spectral resolution of TripleSpec is R = 2500 - 2700, 
corresponding to a Gaussian dispersion ~50 km s . As 
the lowest optical cr* is 67 km s" 1 (for NGC 7457), this 
spectral resolution is suitable for our study. 

For sky subtraction, in particular for the strong OH sky 
emission lines, we also observed blank sky, offset by sev- 
eral arc minutes from each galaxy, since the size of each 
galaxy is larger than the slit length and fills the entire 
slit. We divide the total exposure time into segments of 
200 second exposures to avoid saturation in the if-band. 
The total on-source exposure time ranges from 600 to 
1000 seconds depending on the magnitude of individual 
galaxies (see Table 1). We observed several A0V stars 
each night to correct for telluric lines. We also observed 
11 K- and M-type giant stars as velocity templates for 
the stellar velocity dispersion determination. 

2.2. Data Reduction 

We performed standard data reduction, i.e., bias sub- 
traction, flat-fielding and wavelength calibration using a 
series of IRAF scripts, then extracted one-dimensional 
spectra using various extraction windows. For telluric 
absorption correction, we constructed a telluric line tem- 
plate for each observing night, based on the spectra of 
A0V stars observed during the night. For each A0V star, 
we fitted their Brackett lines with double Gaussians and 
normalized the spectra by its continuum. Dividing the 
observed A0V star spectrum by this model spectrum pro- 
vides a telluric template. We combined all telluric tem- 
plates to construct a mean template for a given night. 
Then, we used the template to correct the galaxy spec- 
tra for telluric absorption lines. 

To investigate the effect of galaxy rotation on stellar 
velocity dispersion measurements, we extracted spatially 
resolved spectra from a number of small extraction win- 
dows (4-16 pixels) along the galaxy's major axis, which 
were allowed to overlap with each other and to slightly in- 
crease at larger radii for obtaining better signal-to-noise. 
The resolution for spatial binning depends on the dis- 
tance to each galaxy. Typically ^10 spectra were ex- 
tracted within a fraction of the R e (see Section 3.5 for 
details) . 

A series of single-aperture spectra were also extracted 
using various aperture sizes, in order to investigate the 
aperture effect. Many previous studies used spatially un- 
resolved cr* measurements, which were affected by line 
broadening due to galaxy rotation. This leads to an over- 
estimation of the galaxy's cr*, depending on how much 
rotation is included in the extraction aperture. Thus, 
the choice of different aperture sizes can affect the cr* 
measurement as presented in § 3.4. In contrast, we can 
correct the cr* measurement for the rotational broaden- 
ing using our spatially resolved spectra. Details on this 
correction are presented in § 3.4 and § 3.5. 

3. ANALYSIS 

3.1. Stellar Velocity Dispersion Measurements 

We measured the stellar velocity dispersion of 31 galax- 
ies in the sample using the stellar lines in the 1.57 - 
1.72 (im range, i.e, CO(4-l) 1.58/mi Si I 1.59^m CO(5- 
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Sample Selection and Observing Log 
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Note. — Col. (1): NGC galaxy catalogue name. Col. (2) Right Ascension. Col. (3): Declination. Col. (4): morphological types. Col. (5): 
distance. Col. (6): spatial scale. Col. (7)-(8): effective radius and reference. Col. (9)-(10): black hole mass and their reference. Col. (11): 
observation date. Co l. (12): total exp o sure ti me. Col. (13): averag e sign al-to-noise r atio w ithin +5 pixel apertu r e. Col. (14): Position Angle. 
Refere nces. — fl) IMarconi Hunt! J20"03T); (2) ICrahaml 12003) : (3) ISani et al.l 12011T): (4) IGultekin et~aH H2009af); (5 ) IMcConnell eTa! 
1 2011]) ; (6) Vcrolmc ct al. (20021) (7) ISchulze L Gcbhar dl II201H); (8) [Bower et al.l H2001I); (9) ILodato fc Bertinl H2003T); (10) ISarzi et a" 
(2001); (11) IDevereux et al. | H2003I); (12) Emscllcm. Dcionghc fc Bacon H999I); (13) IBarth et al.l (120011); (14) Ivan d en Bosch fc dc Zeeuw 
(2010); (15) Giiltckin ct al. (2009b); (16) Hcrrnstcin ct al. (2005); (17) Fcrrarcsc. Ford & Jaffc (1996); (18) Crctton & van den Bosch 
(1999); (19) Walsh, Barth & Sarzi (2010); (20) Gcbhardt ct al. (2011); (21) Shcn & Gcbhardt (2010); (22) Fcrrarcsc & Ford (1999); (23) 
Ivan der Marel fc van den Boschl l |199SI) ; (24) ITremaine et al.1 fl200a) ; (25) IMcConnell fc Ma I ( 120131) 



2) 1.6^m CO(6-3) 1.62L(m CO(7-4) 1.64/im CO(8-5) 
1.66/im, and Mg I 1.7lLjm (see Figure 1). Using the 
Gauss-Hermite Pixel Fitting software (|van der Marell 
[T99llWoo et~all[2 004. 20051 I2006I) . we performed X 2 min- 
imization in fitting the galaxy spectra directly in pixel 
space to stellar template spectra broadened by a Gaus- 
sian kernel with velocity widths ranging from 50 to 350 
km s -1 . The continua of the spectra of the template 
stars are fitted with low-order (2-3) polynomials while 
the Fell emission line at f .65^m, bad pixels and residu- 
als from sky line subtraction were masked out before the 
fitting. 

3.2. Template Mismatch 

Since the ct* measurement is affected by the choice of 
template star, it is necessary to quantify the uncertainty 
due to the template mismatch. Using 11 velocity tem- 
plate stars of various spectral types, namely, K0 III, Kl 



III, K2 III, K5III, two MO III, Ml III, two M2 III, M3 
III, and M5 III, which were observed with the same in- 
strumental setup during our observing runs, we measured 
and compared ct* for individual galaxies in the sample, in 
order to investigate the variation in the ct* measurement 
caused by template mismatch. Then, we accounted for 
template mismatch in the determination of ct* by aver- 
aging ct* measurements using various template stars. 

To compare the overall spectral shapes, we present the 
spectra of the individual template stars, after broadening 
them with a Gaussian velocity (red thick lines) in Figure 
[U The observed spectrum of NGC 1023 is overplottcd 
(black lines) to demonstrate the template mismatch. In 
the stellar spectra, the line strength of the CO absorption 
lines increases toward later-type stars (from upper panels 
to lower panels). This trend is in particular clearly shown 
for the CO (6-3) line at 1.62 /im and also for the Mgl line 
at 1.71 fim. On the other hand, the Sil line strength 
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Fig. 1. — Normalized spectra of template stars. The template 
star spectra (thick red solid line) are compared with the spectra of 
NGC 1023 (thin black solid line). The template star is broadened 
with a Gaussian velocity kernel. Individual template stars show 
different line strengths, particularly for the CO absorption line. 
The K type star templates provide a poor fit to the galaxy spec- 
trum. In the top panel, we marked several individual stellar lines 
with black tick masks (from left, CO(4-l), Sil, CO(5-2), CO(6-3), 
CO (7-4), CO (8-5) and Mgl). 



shows no strong variation with spectral type. The com- 
parison in Fig. 2 clearly shows that spectra of K type 
stars provide a poor match to the observed galaxy spec- 
trum in this wavelength range, while spectra of M type 
stars can fit the observed galaxy spectrum reasonably 
well. 

In Figure^ we compare the multiple measurements of 
<7* of NGC 1023, using each template star for the fitting. 
As expected from Figure 1, the tr* measured from M- 




K0 Kl K2 K5 M0[ M0 2 Ml M2, M2-, M3 M5 

Template spectral type 

Fig. 2. — Comparison of the measured stellar velocity dispersion 
of NGC 1023 using different template stars. Filled circles denote 
M0 III, Ml III, M2 III and M3 III type stars, which were used 
for calculating the mean stellar velocity dispersion. Open circles 
denote M5 III and K type template stars. The mean <r» is given by 
the red diamond and the standard deviation of the measurements 
is included in the uncertainty denoted by the blue solid error bar. 
The blue dashed line indicates the range of template stars used for 
calculating the mean stellar velocity dispersion. 



type stars shows small variation, while er* measured from 
K-type stars exhibits a larger scatter, suggesting that M- 
type stars provide a fair representation of the luminosity- 
weighted stellar population in the _ff-band. Therefore, we 
excluded the measurements from K-type stars and calcu- 
lated the mean based on the 6 M-type stars. Since the 
M5 III star shows slightly different line shapes compared 
to other M-type stars (see Figure 1), we also excluded 
the measurement based on the M5 III template. After 
calculating the standard deviation of the measurements 
from 6 M-type stars as the uncertainty of template mis- 
match, we added the uncertainty of template mismatch 
to the mean measurement errors from 6 M-type stars in 
quadrature, in order to determine the uncertainty of <7*. 
For example, the red diamond in Figure [2] indicates the 
mean er*, derived from 6 template stars and its uncer- 
tainty. In Figure [31 we present the normalized observed 
spectrum (black solid line) of each galaxy in the sample, 
overplotted with the best-fit model (red solid line). 

3.3. Spatially Resolved Stellar Velocity Dispersions 

By extracting spectra over several bins along the major 
axis, we obtained spatially resolved kinematics. In Fig- 
ures H] and [5] we show the radial profiles of line-of-sight 
velocities (upper panel) and velocity dispersions (lower 
panel) for each galaxy. We used the line-of-sight veloc- 
ity of the galaxy center as a reference and normalized 
all velocities with respect to the central value. For most 
galaxies we extracted 9-13 spectra along the slit (in the 
direction of the major axis) out to ±7" from the center. 
This is smaller than the slit size (±15") since we were not 
able to use outer pixels due to much lower S/N than the 
central bins to measure and the ABBA dither pattern 
along the slit. 

Among the sample galaxies, we find a clear rotation 
component for 25 out of 31 objects. The amplitude of 
the projected rotation velocity ranges from ^40 km s _1 
to over 200 km s _1 while six galaxies, namely NGC 1068, 
NGC 3608, NGC 4261, NGC 4374, NGC 4486 and NGC 
6251 show a weak or no rotation component. For the 
galaxies with a significant rotation component, we expect 
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Fig. 3. — Normalized spectra of the 31 galaxies and their best-fit 
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well. Residuals of OH sky emission lines (sharp fcautures in Fig. 1) 
and the AGN Fe II emission line (e.g., in NGC 1068) were masked 
out before fitting. 

line broadening due to the rotation, leading to overesti- 
mation of <7* , if a large single aperture is used for extrac- 
tion. While the magnitude of this effect depends on the 
details of radial profiles of rotation and velocity disper- 
sion of each galaxy, it will lead to a systematic bias if not 
taken into account. We will account for the rotation com- 
ponent in the stellar velocity dispersion measurement in 
the next subsection. 

In the lower panels of Figure 0] and Figure [5j we show 
the stellar velocity dispersion profiles along the major 
axis. While stellar velocity dispersion decreases from 
the center to the outer regions for most galaxies, sev- 
eral galaxies, e.g., NGC 1068, NGC 4261, NGC 4374, 
NGC 4596 and NGC 7052, do not show such a decreas- 



ing trend of cr*, but rather show irregular shapes; fiat, 
increasing or asymmetric trends as similarly reported 
by prev ious studies based on optical kinematics stud- 
ies (e.g. IDresslerl 119841: IBender. Saglia fe Gerhard1ll994t 
[Kentj|1990l : IPinknev et al.ll2003f) . 



3.4. Aperture Size Effect 

We investigate the effect of using different aperture 
sizes on the measured stellar velocity dispersion, by di- 
rectly measuring cr* from apertures of increasing size. In 
Figure|5]we present the cr* measurements as a function of 
aperture size, after normalizing them to the cr* measured 
from the smallest aperture (4-10"). We find three differ- 
ent trends (increasing, flat and decreasing) of cr* with 
increasing aperture size. The 8 galaxies shown in the 
upper panels exhibit an increase of larger aper- 

ture size is used, while for 13 galaxies shown in the lower 
panels, cr* decreases with increasing aperture size. These 
galaxies show variation of cr* up to 20% as aperture size 
changes. In contrast, 10 galaxies (middle panels) do not 
show clear change of cr* as a function of aperture size. 

Thus, when measuring cr* from a large aperture it is 
possible to either overestimate or underestimate cr* . The 
magnitude and direction of this bias depends on two fac- 
tors: (1) the overestimation caused by rotational line 
broadening and (2) the natural decrease of func- 
tion of radius. The galaxies in the upper panels in Fig- 
ure [5] are dominated by the first effect. They show rel- 
atively strong galaxy rotation and only a mild decrease 
in their velocity dispersion profile, as shown in Figures 
[4] and [5l leading to a net increase in cr* with increasing 
aperture size. For these galaxies, when the extraction 
aperture covers outer parts of the galaxy, where the ro- 
tation curves flatten, the aperture effect on cr* also flat- 
tens. This is clearly seen for example in the case of NGC 
3384, which shows a flattening of cr* beyond the third bin, 
corresponding to the flattening of the rotation curve in 
Figure [4] NGC 4742 even shows a decrease in cr* beyond 
the radius where the rotation curves becomes flat. 

The decreasing cr* trend for the galaxies in the lower 
panels of Figure [5] can be explained in a similar way. 
For these galaxies the decrease in cr* profile is domi- 
nating over line broadening due to rotation. For exam- 
ple, NGC 4486 shows the largest variation in cr* as a 
function of aperture size since it has no significant ro- 
tation component while the velocity dispersion profile is 
strongly decreasing toward larger radii. Similarly, NGC 
3608 and NGC 4649 also show strong decrease, which is 
dominated by the strong decrease in cr* with radius. In 
contrast, there are galaxies, e.g., NGC 821, where the 
effect from the rotation curve and the decreasing cr* pro- 
file are roughly of the same order, leading to a small net 
variation of cr* measured from different apertures. 



3.5. Correction of the Galaxy Rotation Effect 

As discussed above, cr* measured within a certain aper- 
ture will be susceptible to line broadening by galaxy ro- 
tation. In contrast, the spatially resolved stellar velocity 
dispersions, represented in Figure 2] and Figure [SJ do 
not suffer rotational broadening. Thus, we can use these 
measurements to compute a rotation-corrected cr*. We 
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Fig. 5. — Same as Figure|4]for the rest of the sample. 
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compute luminosity-weighted er* within a radius R: 



f R R v*(r)I(r)dr 
f« R I(r)dr 



(1) 



where I(r) is the surface brightness profile of the galaxy 
and R is the outer radius within which we compute cr*. 
Using the spectral images, we measure the surface bright- 
ness profile of each galaxy by fitting with two Gaus- 
sian models, and use this fit to compute the luminosity- 
weight for the <7* measured at each radius. We chose 
an outer radius R for each galaxy based on the R e (see 
Table 1). For 20 galaxies, we were able to measure spa- 
tially resolved c* only at the central parts, due to the 
limited spatial coverage and/or lower S/N at the outer 
part. Thus, we chose 1/8 of R c as an outer radius in 
Eq. 1. For the other 11 galaxies, we measured er* over a 



larger fraction of R c (1/4 to unity) as an outer radius and 
corrected for the rotation component as listed in Table 
2. 

In Figure [7] we illustrate the effect of the correction for 
galaxy rotation. Here, an is the luminosity-weighted cr» 
within R as computed from Equation [T] while an(uncor) 
is measured from a single aperture with an aperture size 
of R. As expected, most galaxies show a decrease in ve- 
locity dispersion when accounting for the rotation com- 
ponent while for galaxies without strong rotation com- 
ponent, the correction is marginal. Including 6 galax- 
ies that show no rotation, the average correction is 6%, 
while the correction for individual galaxies can be up to 
^20%. The magnitude of the rotation correction tends 
to be smaller for more massive galaxies. NGC 7052 
with the highest cr„ in the sample seems to be an outlier 
from this trend since it has relatively large rotation while 
cr, mildly decreases within R c /8. In summary, we find 
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Fig. 7. — Ratio of velocity dispersions determined with and with- 
out rotation correction as a function of velocity dispersion. Filled 
and open circles represent elliptical and lenticular galaxies, respec- 
tively, while stars denote late-type galaxies. Six galaxies without 
a clear rotation component are marked with open squares. 

that stellar velocity dispersions measured from single- 
aperture spectra can be bias ed by up to ~ 20%. T his is 
consistent w i th the results of lBennert et al.l (|2011al ) and 
lHarris et all (|2012t ). 

A possible drawback for the comparison with previous 
studies, that usually report cr* measured within R e , is the 
limited spatial coverage in our work, restricting our mea- 
surements to Re/8. To investigate the effect on the cr* 
measurement, we tested two extreme cases. First, we as- 
sumed a constant cr* from R e /8 to R c , equal to the value 
at Re/8. Second, we extrapolated the decreasing stel- 
lar velocity dispersion profile out to R c . For both cases 
we computed the luminosity-weighted cr* within R e via 
Equation Q] We found that cr* values decrease by only 
a few per cent by increasing the outer radius from R e /8 
to R e in both cases, due to the much lower luminosity 
weight at outer radii. Thus, our cr* measurements within 
Re/8 will closely resemble the value that would be mea- 
sured at R c . 

4. DISCUSSION 
4.1. Optical vs. near-IR Velocity Dispersions 

In Figure 8 we compare the stellar velocity dispersions 
measured using our near-IR spectra (am) with the lit- 
erature values measured from optical spectra (a op t). For 
this comparison, we col lected a op t measurements from 
iMcConnell fc Mai (|2013l ). who listed their own measure- 
ments as well as previously measured values from the 
literature (see Table 2). Although these optical measure- 
ments were based on spatially resolved stellar kinematics 
and the quoted values were luminosity-weighted or aver- 
aged velocity dispersions within the R e for most galaxies, 
these values were not homogeneously measured due to 
the various data quality and the measurement methods. 
In addition, some of these values in the original works 
were measured with a smaller aperture size than R e or 
the aperture size was not clearly stated for many cases 
in the original references. 

Thus, we decide to compare both //-band cr* measure- 
ments with/ without rotation correction to the optical cr* 
measurements. As shown in Figure [8j we find no sig- 



TABLE 2 

Near-IR and optical stellar velocity dispersions 



Name ctir cr op t Ref. 

(km s" 1 ) (km s" 1 ) (R c ) (km s" 1 ) (km s" 1 ) 
(1) (2) (3) (4) (5) (6) (7) 



Note. — Col. (1): Object name. Col. (2): //-band stellar velocity 
dispersion measured using a single aperture size R without rotation 
correction. Col. (3): luminosity- weighted //-band a* within R. Col. 

(4) : Aperture radius R used for //-band cr, in units of R c . Col. (5): 
//-band velocity dispersion measured using a single aperture within 
±7". Col. (6): optical stellar velocity dispersions. Col. (7): Reference 

for optical velocity d ispersions. 

Referen ces. - f 1) Ivan der Marel eTaTI 11991 1 ; (2) ICebhardt et al.l 
(I2000al); (31 IPinknev et al.l 1 12003); (4) IBower et al.l l !200lT )7 

(5) INelson fe WhittleT lfl995l) : (6) [Gultckin et al. ( 2009a) ; (7 ) 
IKorme ndv fc Gcbhardt (2001); (8) Kormcn dy et al.l (1998): 

(9) ICebhardt et al.l d200aT : (10) Gcbhard t"et al.l (1200011); (11) 
IGiiltekin et al.l (I2009bl); (12) Ivan der Marel, Binnev & Davies 

111990) ; (13) ICretton fc van den Boschl jjl999l ): ( 14) IBower et al. 

|l998|) ; (15) ICebhardt et alj (120111); (16) IMcConnell etaT 

]201ll); (17) ISmith, Heckman fe Illingworthl (119901 ); (18) 

van den Bosch & van der Marell l |1995l ) 



nificant difference between o~ opt and am. The best fit 
between optical and near-IR measurements is close to a 
one-to-one relation with a scatter of ^0.04 dex (10%) 
when rotation is not corrected for am- The lower panel 
in Fig. 8 shows that the rotation-corrected am is slightly 
smaller than the optical cr* at lower mass range. How- 
ever, the average offset is only 7%, which is not significant 
compared to the measurement uncertainties of stellar ve- 
locity disper sions. 

Note that IMcConnell fc Ma I (|2013f ) included rotation 
in calculating luminosity-weighted cr* by adding rotation 
velocity to velocity dispersion in quadrature (See their 
Eq. 1). Thu s, the slight offset betwee n our am and op- 
tical cr* from IMcConnell fc Ma I (|2013f) may be explained 
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Fig. 8. — Comparison between _ff-band stellar velocity disper- 
sions and optical stellar velocity dispersions from the literature. 
H-band <t» is measured from an extraction aperture, ±7" (top 
panel) or corrected for rotation effect (bottom panel). Elliptical 
and lenticular galaxies are denoted with filled and open circles while 
3 late-type galaxies are indicated by filled stars. The best-fit (solid 
line) is consistent with a one-to-one relation (dotted line). 

by the rotation effect. To test this hypothesis, we derive 
rotation-included velo city dispersions using the same in- 
tegral as adopted by IMcConnell k Ma I (|2013t Eq. 1). 
For these consistently measured velocity dispersions, we 
find that optical and IR velocity dispersions show a one- 
to-one relationship with a slope of 1.00 ± 0.05 and a 0.03 
dex (7%) intrinsic scatter. Thus, we conclude that <r* 
measurements derived from optical and H-band stellar 
lines are consistent. 

While many studies have been devoted to measure c* 
of galaxies using optical spectra, the number of er* studies 
based on near-IR data, either ii-band or if-band spectra 
is growing. However, there are currently only few stud- 
ies that actually compared th e results from both wave- 
length regimes. For example, ISilge k GebhaTdl {2003) 
measured <r* of a sample of 25 elliptical and lenticular 
galaxies using the 2.29/im CO(2-0) band head in the 
.ftf-band spectra. Comparing their IR results to optical 
velocity dispersions from the literature, they concluded 
that IR stellar velocity dispersions can be lower than op- 



tical stellar velocity dispersions, by up to 30% - 40% 
and with a median offset of 11%. The inconsistency be- 
tween optical and near-IR measurements is probably due 
to a sample bias and measurement uncertainties. First, 
their sample mainly consists of SO galaxies and the sys- 
tematic difference between optical and near-IR velocity 
dispersion in their study is mainly caused by SO galaxies 
while their elliptical subsamplc docs not show a differ- 
ence between optical and nea r-IR measurements. Sec- 
ond, ISilge k Gebhardtl (|2003f ) measured cr» using solely 
a single CO band head in the iT-band, which is much 
more susceptible to template mismatch as explained in 
their analysis. 

In contrast t o ISilge k Gebhardtl (|2003ft . 
lRothberg~fc Fischer] (|2010l ) reported no inconsistency 
between optical and near-IR stellar velocity dispersion 
measurements for elliptical galaxies. Using a sample 
of 23 elliptical galaxies and 14 merger remnants, they 
measured <r* from stellar lines in the if-band spectra, 
i.e., CO (2-0), CO (3-1), and CO (4-2) band heads, 
and compare them with velocity dispersion measured 
from the optical Call triplet line, showing that optical 
and near-IR stellar velocity dispersions are virtually 
the same for ellipti cal gala xies. For merger remnants 
iRothberg k Fischer] {2010;) reported a discrepancy 
between optical and near-IR velocity dispersions, pre- 
sumably due to the presence of young stellar population, 
which are obscured at optical wavelengths. However, 
for elliptical galaxies, their results indicate that optical 
and near-IR stellar lines represent the same kinematics 
and a du st effect is neglig i ble. Similarly, a recent 
study by iVanderbeke et al l (|2011[) presented near-IR 
it* measurements also based on the CO band heads for 
a sample of 22 galaxies, consisting of similar numbers 
of ellipticals and lenticulars. Comparing with previous 
optical measurements they reported that optical and 
near-IR c* were consistent for their sample, which is 
consistent with our results. 

In the case of velocity dispersions using ii-band stel- 
lar lines, there has been no systematic comparison with 
optical velocity dispersions. Using various stellar lines 
in the ii-band spectra and carefully accounting for the 
template mismatch problem (see Section 3.2), for the 
first time, we show that optical and ii-band <r* mea- 
surements are consistent for early-type galaxies, indi- 
cating that optical and ii-band stellar lines represent 
the same kinematics and that a dust effect, i.e., ob- 
scuration at optical wavelengths, is negligible. These 
results are consistent with Jf-band stellar kinem atics 
(|Rothberg fc Fischerll2010l: IVanderbeke et alj|201ll ). Our 
results imply that near-IR cr, measurements carried out 
for AGN host galaxies, for which optical measurements 
are more difficult to perform due to the strong AGN 
contribution (e.g., Woo et al. 2010), provide unbiased 
results, compared to optical measurements. 

4.2. The Mbh-c* Relation for Early-type Galaxies 

In this paper we present homogeneously measured c* 
for 31 galaxies with dynamical Mbh measurements. By 
accounting for galaxy rotation and implementing a uni- 
form analysis for measuring velocity dispersions, our er„ 
measurements are slightly different from previous optical 
values. In this section we demonstrate the effect of these 
new velocity dispersions on the Mbh-c* relation by fit- 
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Fig. 9. — The Mbh-c* relation of 31 nearby galaxies, us- 
ing rotation-corrected cr, measured from our H-band spectra 
(red circles, red line) and rotation-included optical cr* from 
IMcCon ncll & Ma (20l3) (gray diamonds, gray line), respectively. 

ting the Mbh-c* relation for 31 galaxies, for which we 
obtained the rotation-corrected cr* . Results on the Mbh- 
cr* relation for the full sample of galaxies with dynamical 
A/bh will be presented in a companion paper (Woo et al. 
2013 in preparation). For this analysis, we used the most 
recent M BH updates from IMcConnell fe Mai (|2013f ). 
We fit the Mbh-c* relation as a single-index power law: 

log(M BH /M Q ) = a + /31og(cr* /200 km s" 1 ) . (2) 

We used the FITEXY method, modified to account for 
intrinsic scatter in the relation ()Tremaine et al.l 120021 : 
iPark et ahll2012| ). to perform the fit as shown in Figure 

m 

By fitt ing the M B h-o* re lation using the optical 
cr* from IMcConnell fe Mai (|2013f ). we obtain a = 
8.27±0.07, P = 3.98±0.40 and an intrinsic scatter of 
0.38±0.05 dex. Our sample has a large overlap with the 
sample used by iTremaine et al.l (|2002f ) . and indeed we 
obtain consistent results for the M B h-c* relation. The 
slight difference arises mainly from a few different galax- 
ies in both samples and from updated Mbh determi- 
nations (e.g. ISchulze fc Gebhardtll20ll . By fitting the 
Mbh-c* relation using our rotation-corrected near-IR cr*, 
we find a = 8.35±0.08, /? = 3.78±0.43 and an intrinsic 
scatter of 0.38±0.05 dex. This relation is slightly shal- 
lower than that derived from optical cr*, but consistent 
within the uncertainties. 

Note that the previous studi es on the Mbh-c* relation 
bv lGiiltekin eFatl (I2009al) and IMcConnell fc Mai (l2013h 
explicitly included rotation in calculating luminosity- 
weighted cr*. In this case, we expect systematic effect on 
the measured cr* due to the random orientation of stellar 
disk with respect to the line-of-sight. To investigate this 
effect, we calculated cr* by adding velo city to velocity dis- 
persio n in quadrature using Eq. 1 in IMcConnell fc Mai 
(|2013l) . The result shows that rotation- included cr* is 



slightly larger than rotation-corrected cr*, particularly 
at low mass scale, by 0.02 dex (~5%) on average with 
a 0.02 (^5%) scatter. Consequently, when we replace 
rotation-corrected cr* with rotation-included cr* in fitting 
the the Mbh-c* relation, the slope slightly increases from 
3.78 ±0.43 to 3.97 ±0.49 as the cr* values increase prefer- 
entially at low mass scale, while intrinsic scatter remains 
the same. 

Although we expect that the rotation effect will sys- 
tematically affect the M B h-c* relation, we do not clearly 
detect the improvement of the Mbh-c* relation by using 
rotation-corrected cr*, presumably due to two reasons. 
First, rotation effect on the luminosity-weighted velocity 
dispersion may not be significant as the integrated veloc- 
ity dispersions are dominated by the inner part, where 
the rotation velocity is relatively small. In the case of 
the rotation-included velocity dispersion, the luminosity 
weight of the inner part is more dominant since velocity 
dispersions are integrat ed in quadrature (see Eq. 1 in 
IMcConnell fc Mai 120131) . Secondly, since our sample is 
mainly composed of early- type galaxies, rotation effect 
is relatively weak compared to late- type galaxies. For 
late-type galaxies with a low cr*, rotation effect can be 
significant, hence, it would be essential to correct for, in 
order to properly derive the Mbh-c, relation. 

For massive BHs, the sphere of influence of BH can 
be large enough to change the effective cr* measurements 
since the velocity dispersion at the center increases due 
to the presence of a BH. Thus, by excluding the sphere 
of influence of BH in calculating the luminosity- weighted 
cr*, the effective velocity dispersions will be decreased. 
Since these corrections can be done only for massive 
galaxies with a resolved sphere of influence, the slope 
of the Mbh-c* relation will increase due to the prefer- 
ential decrease of cr* at high mass end. For example, 
by excluding the sphere of influence of BH in deriving 
the effective cr* wi t hin R e for 12 most massive galaxies, 
IMcConnell fc Ma I (|2013l) showed that the slope of the 
M B H-cr* relation increased from 5.48±0.30 to 5.64±0.32. 
We performed a similar analysis using our data although 
only two galaxies, NGC 4486 and NGC 4649, are among 
those 12 galaxies with a resolved sphere of influence. By 
excluding the sphere of influence of BH, the luminosity- 
weighted cr* decreases from 346 ± 12 to 327 ± 11 for NGC 
4486, and from 346±12 to 327±11 for NGC 4649. Based 
on these two updated cr* , the slope of the Mbh-c* rela- 
tion slightly increases from 3.78±0.43 to 3.79±0.45, how- 
ever two slopes are consistent within the uncertainties. 
Using only 2 galaxies, it is not clear whether excluding 
or including the sphere of influence in determining the 
effective cr* improves the Mbh-c* relation. 

Compared to the Mbh-ct * relation recently presented 
bv lMcConnell fc Ma I (|2013h . we find a significantly shal- 
lower slope. They report a slope of 5.64 ± 0.32, using a 
much larger galaxy sample, which includes in particular 
more galaxies at higher and lower masses. We will in- 
vestigate the implications of our results on the Mbh-ct* 
relation in detail in a companion paper (Woo et al. 2013 
in preparation). 

5. SUMMARY 

We observed a sample of 31 nearby galaxies with 
Triplespec, a near-IR long-slit spectrograph at the Palo- 
mar 5-m telescope in order to homogeneously measure 
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velocity dispersions from the H-band stellar lines. The 
galaxies in the sample cover a wide range in (67 km 
s _1 < a* < 385 km s _1 ) and their dynamical central BH 
masses are also available. To account for template mis- 
match, we used 11 giant stars with spectral type ranging 
from KO to M5 as velocity templates, and found that M 
giants generate the most reliable fits and velocity disper- 
sion measurements. 

By measuring velocity and velocity dispersion as a 
function of radius along the major axis of each galaxy, 
we determined the rotation curve and velocity dispersion 
profile. Using these spatially resolved velocity dispersion 
measurements, we calculated the luminosity- weighted 
stellar velocity dispersions within the R c of each galaxy. 
For 25 out of 31 galaxies in the sample, we found a clear 
rotation component, indicating that stellar velocity dis- 
persions can be significantly overestimated due to the 
rotational broadening if a large single aperture is used to 
extract spectra. Compared to rotation-corrected velocity 
dispersions, velocity dispersions measured from single- 
aperture spectra showed systematically larger values by 
up to -20%. 

We compared velocity dispersions measured from H- 
band stellar lines with those measured from optical lines 
and found no systematic difference, suggesting that op- 
tical and .ff-band stellar lines represent the same kine- 
matics and that dust effect is negligible for early-type 
galaxies. Our results confirm that optical and near-IR 
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stellar lines can be interchangeably used to measure stel- 
lar kinematics and near-IR cr, measured for AGN host 
galaxies can be directly compared to optical ct* of quies- 
cent galaxies. 

Using the rotation-corrected ct» measurements based 
on the spatially- resolved H-band spectra of 31 nearby 
galaxies, we derived the Mbh-c* relation to investi- 
gate rotation effect. The slope of the Mbh-c* relation 
is slightly shallower than that based on the rotation- 
included optical or near-IR ct* measurements. Although 
rotation effect is not dramatically strong for early-type 
galaxies, it is potentially important to correct for, par- 
ticularly for low mass, late-type galaxies with a strong 
rotation component, in order to properly determine the 
Mbh-c* relation and its intrinsic scatter. A future study 
based on spatially resolved spectra for late-type galaxies 
is required to fully quantify rotation effect on the Mbh-c* 
relation. 
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